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Abstract. After more than 15 years of quiescence the Be/X-ray binary V0332-r53 underwent a giant outburst in December 
2004. We have investigated the timing properties of the source in correlation with its spectral states as defined by different 
positions in the colour-colour diagram. We have used RXTE and INTEGRAL light curves to obtain colour-colour diagrams, 
power spectra and phase-lag spectra. The power spectra were fitted with a multi-Lorentzian function. The source shows two 
distinct branches in the colour-colour diagram that resemble those of the Z sources. The hard branch (similar to the horizontal 
branch of Z sources) is characterised by a low-amplitude change of the hard colour compared to the change in the soft colour. In 
the soft branch (analogue to the normal branch) the amplitude of variability of the hard colour is about three times larger than that 
of the soft colour. As the count rate decreases the source moves up gradually through the soft to the hard branch. The aperiodic 
variability (excluding the pulse noise) consists of band-limited noise (represented by three broad Lorentzian components) and 
two QPOs at 0.05 Hz and 0.22 Hz. The strength of the lower frequency QPO increases as the source approaches the hard branch 
(similar to HBOs in Z sources). The higher frequency QPO reaches maximum significance when the source is in the middle 
of the branch (like NBOs). We have performed the first measurements of phase lags in the band limited noise below 8 Hz 
in an accreting X-ray pulsar and found that soft lags dominate at high frequencies. Above the pulse frequency (0.23 Hz), the 
amplitude of the lag increases as the X-ray flux increases. The Z topology appears to be a signature of the neutron star binaries 
as it is present in all types of neutron-star binaries (Z, atoll and, as we show here, in accreting pulsars as well). However, the 
motion along this track, the time scales through the different branches of the diagram and the aperiodic variability associated 
with portions of the Z track differ for each subclass of neutron-star binary. 

Key words, stars: individual: V0332-t54, BQ Cam - X-rays: binaries - stars: neutron - stars: binaries close -stars: emission 
line. Be 


1. Introduction 


V0332-I-53 is a hard X-ray transient, accreting X-ray pulsar 
and Be/X-ray binary that spends most of its life in a quies¬ 
cent X-ray state. This quiescent state is occasionally inter¬ 
rupted by sudden increases of the X-ray flux in which the 
source reaches Eddington luminosities. V0332H-53 was discov¬ 
ered during one of these large outbursts in 1973 by the Vela 
5B satellite (iTerrell & Priedhorskvll984HWhitlockll989l) . The 
outburst lasted for about 100 days and reached 1.4 Crab, or 
1.5 X 103« erg s ' in the 3-12 keV energy range, assuming a 
distance of 7 k pc jNegueruela et al.lll999l) . V0332-I-53 reap¬ 
peared in 1983 ( iTanakal 19831) in the form of three small out¬ 
bursts. EXOSAT observations of this activity period resulted 
in the discovery of X-ray pulsations with Pspin = 4.4 s and 
the determinati on of the orbital parameters, Porh - 34.25 d 
and e = 0.31 dStella et alJll985l) . Ginga detected V0332-I-53 
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again in 1989 with a peak flux of 0.4 Crab (1-20 keV). The 
analysis of this new outburst with improved technology al¬ 
lowed the discov ery of a cyclotron reso nant scattering fea¬ 
ture at 28.5 keV ilMakishima et al.lll990h and QPOs at 0.051 
Hz JTakeshima et al.lll994h . The last outburst took place in 
2004 wit h a peak intensity in the 1.3-12.1 keV range of ~ 
1.1 Crab iSwank et al.l2004llRemillardl2004 . INTEGRAL ob¬ 
servations provided the first broad- band spectrum (5-100 keV ) 
and detected three cyclotron lines ilKrevkenhohm et ^l200.5[) . 
Using RXTE obse rvations obtained during the 2004 outburst, 
|Zhang_e^ayJ200^refined the orbital parameters of the system 
and lOu eTalT(l20(^ reported the discovery of a new QPO at 
0.22 Hz. 


The optical counterpart to V0332-I-53 is an 08-9Ve star at 
a distance of ~ 7 kpc, show ing Her in emission and strong and 


variable infrared emission dRernacca et al .111 984 

; ICorbet et all 

ll986HCoe et al.ll987HNegueruela et alJl999). 
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Fig. 1. The outburst of V0332+53 as seen by ASM RXTE (1 
Crab ~ 75 c/s). Grey and black lines mark the time of the RXTE 
and INTEGRAL observations, respectively. 

In this work we have analysed data from the JEM-X 
INTEGRAL and the PCA RXTE instruments during the last 
large 2004 outburst. Our aim is to perform a correlated tim¬ 
ing an spectral analysis of the X-ray emission as the outburst 
decayed. 

2. Observations and data analysis 

We have used data from the space missions Rossi X-ray 
timing Explorer (RXTE) and the INTEmational Gamma-ray 
Astrophysics Laboratory (INTEGRAL). The RXTE data cor¬ 
respond to public Target of Opportunity (TOO) observations 
performed between December 28, 2004 and March 17, 2005 
(MID 53368.24-53447.06). Given the high count rate and the 
fact that the PCU number 2 was on all the time we used data 
from this detector only. Data reduction was performed with 
ETOOLS v5.3.L The total RXTE on-source time amounted to 
~ 113 ks. 

The INTEGRAL data were retrieved from the public 
archive and correspond to the Target of Opportunity (TOO) 
observations performed between January 5, 2005 and Eebruary 
21, 2005 (MJD 53376.27-53422.05). We analysed JEM-X data 
for a total observing time of ~ 338 ks, divided in 105 point¬ 
ings (Science Windows) in revolutions 273-274, 278, 284- 
288. Data reduction was carried out with the standard Off-line 
Scientific Analysis (OSA) software version 5.0, available from 
the INTEGRAL Science Data Centre (ISDC)'. 

The timing analysis was done with the XRONOS software 
package v5.21 jStella & Angelinill992h . 

3. Results 

3.1. The outburst 

The X-ray outburst began in 2004 November (MJD 53330) 
and reached maximum flux about one month later (~ MJD 
53368). While RXTE observations coincided with the peak of 
the outburst, INTEGRAL observations started in 2005 January 
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5 (MJD 53376), i.e., a few days after the peak. Given that the 
observations only covered the decay of the outburst and in or¬ 
der to have a clear picture of the outburst profile we show the 
RXTE ASM light curve in Eig. [Q The total duration of the 
outburst was ~ 100 days, of which the INTEGRAL observa¬ 
tions covered about 46 days and the RXTE observations the 
last 80 days (from maximum flux to quiescence). The peak 
flux, as measured by the ASM (1.3-12.1 keV) was ~L1 Crab 
lRemillardi2004l) . 

The outburst profile is very symmetric, although the rise is 
somewhat faster than the decay. The difference appeared, how¬ 
ever, at the end of the decay, with the profile displaying a longer 
tail at the base of the outburst. While a decrease in flux of about 
30% is achieved in ~ 5 days at the onset of the decay, such a 
drop requires 12 days at the end of the decay. Surprisingly, this 
behaviour contrasts with that observed in A 0535-1-262 during 
its June 2005 Type II outburst, where the longe r tail was seen 
at the beginning of the outburst jCoe et al .1200511 . 

At the beginning of the RXTE observations the 4-30 keV 
PCU2 count rate was ~ 3000 c/s, that is, ~L5 Crab. The in¬ 
tensity of the last observation reported here was 37 c/s, which 
corresponds to ~0.02 Crab. Overall the decrease in flux was a 

factor ~ 60, from 5.3 x 10“^ erg cm“^ s“^ to 0.09 x 10“^ erg 

-2 -1 
cm s . 

In order to study the evolution of the timing and spectral 
properties of V0332-I-53 throughout the decay of the outburst 
we divided the INTEGRAL and RXTE light curves into several 
regions. The start time of these regions have been marked with 
black (INTEGRAL) and grey (RXTE) lines in Eig.H Tabled 
gives the start and exposure time for each region. 

3.2. Colour-colour diagram 

Background-subtracted light curves corresponding to the en¬ 
ergy ranges ci = 4 - 7.5 keV, C 2 - 7.5 - 10 keV, C 3 = 10 - 15 
keV and C 4 = 15 - 30 keV were used to define the soft colour 
(SC) as the ratio ca/ci and the hard colour (HC) as the ratio 
between C4/C3. The colour-colour diagram (CD) of V0332H-53 
was then constructed by plotting the hard colour as function of 
the soft colour (Eig. Ell. 

At the peak of the outburst the source was in a soft state. 
As the count rate decreased the source moved up and right 
in the CD, i.e. it became harder. At the end of the outburst 
(INTEGRAL regions H and I and RXTE regions E and E) the 
source moved to the left, that is, toward lower values of SC 
with little variation in the HC. Two spectral states or branches 
can be distinguished in the CD: a soft state that corresponds 
to higher count rates and a hard state that includes points with 
lower count rates. In the soft branch, the main parameter driv¬ 
ing the spectral variability is the hard colour: while the HC 
changes by about 45%, the SC changes by about 15%. In the 
hard branch, the variability is more pronounced in the SC. 

The pattern traced out in the CD by V0332-H53 rese mbles 
that of the low-mass Z sources (see e. g ivan der Klisll995l) . with 
the soft branch being the analogue to the normal branch and 
the hard branch being the counterpart of the horizontal branch. 
The flaring branch would be missing in V0332-I-53. As in Z 
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Fig. 2. JEM-X and PCA observations of the outburst. Flux evolution (top) and colour-colour diagrams (bottom). 1 Crab corre¬ 
sponds to ~ 2000 PCU2 c/s and ~130 JEM-X c/s in the energy range 4-30 keV. The X-ray colours are defined as SC=7.5-10 
keV/4-7.5 keV and HC=15-30 keV/10-15 keV, where the energy ranges represent the background subtracted count rates in the 
corresponding band. 


sources, the hardest spectrum corresponds to a lower count rate 
state. Also, the source moves gradually through the different 
branches without jumps. The source took about 80 days to go 
from the softest point of the soft branch to the hardest point of 
the hard branch. 

3.3. Power spectral analysis 

In order to investigate the aperiodic variability of V0332-I-53 
in relation with its spectral state we divided the CD into six 
regions and calculated power spectra for each region. The 
power spectra were obtained by dividing the 0.125-s resolu¬ 
tion PCA light curves into 256-s segments and calculating the 
Fast Fourier Transform (FFT) of each segment. We obtained 
a mean power spectrum for each region by taking the average 
over the segments included in the region. See Table[2to find the 
average values of the colours and intensity of each region. The 
power was normalised such that the integral gives the squared 
nns fractional variability. Given the superior quality, in terms 
of the signal-to-noise ratio, and the longer outburst coverage of 
the RXTE observations we only used RXTE data in the timing 
analysis. 

Figure 0 shows the power spectra associated with dif¬ 
ferent positions in the CD. In order to fit the power spec¬ 
tra we followed the approach commonly employed recently 
for low-mass neutron-star and black-ho l e systems of using 
Lorentzian functions only (lNowakll200fA iRelloni et al.ll20o3 


Table 1. Results of the timing and colour analysis 


Reg. 

Start time 

MJD 

Exposure 
time (ks) 

Count 

rate“ 

SC 

HC 

Flux"* 

INTEGRAL data 

A 

53376.27 

32.842 

229 

0.76 

0.60 

4.3 

B 

53379.06 

42.656 

216 

0.77 

0.63 

4.2 

C 

53380.22 

10.504 

194 

0.79 

0.64 

3.9 

D 

53393.64 

52.193 

104 

0.80 

0.70 

2.1 

E 

53410.14 

138.368 

53 

0.83 

0.84 

1.1 

F 

53413.07 

13.316 

46 

0.85 

0.86 

0.9 

G 

53416.45 

15.403 

33 

0.85 

0.88 

0.6 

H 

53419.05 

14.725 

23 

0.84 

0.96 

0.4 

I 

53422.05 

17.900 

14 

0.74 

0.95 

0.3 

RXTE data 

A 

53368.24 

6.688 

2974 

0.52 

0.21 

5.3 

B 

53374.03 

5.520 

2566 

0.54 

0.23 

4.8 

C 

53384.36 

14.928 

1742 

0.55 

0.24 

3.4 

D 

53386.86 

41.984 

1464 

0.56 

0.25 

2.9 

E 

53413.06 

38.608 

406 

0.60 

0.31 

0.9 

F 

53436.30 

5.120 

37 

0.48 

0.40 

0.09 


a: background subtracted 4-30 keV 
b\ 10“* erg cm“^ 


Ivan Straaten et al .120021) . Note that with the FrequencyxPower 
representation used in Fig.|3 the central frequency, Vc, of the 
Lorentzian function and the frequency at which the Lorentzian 
reaches maximum power, Vmax, do not coincide but are related 
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Frequency (Hz) 


Fig. 3. Evolution of the RXTE power spectra during the decay of the outburst. The various Lorentzian components have been 
indicated: Lvlfn (dot), Lhfn (dash-dot-dot), Llfqpo (dash), Lhfqpo (dash-dot) and Llfn (dash-dash-dot). 


by Vmax = (Vc + where A is the Lorentzian half-width at 

half-maximum. 

The study of the aperiodic variability in V0332-F53 is ham¬ 
pered by the presence of very narrow peaks in the power spec¬ 
tra. These spikes correspond to the spin period of the system. 
V0332-F53 is an X-ray pulsar with Pgpm = 4.37 s, which corre¬ 
sponds to a frequency of 0.228 Hz. This peak and up to four 
harmonics can be seen in the power spectra (see Eig. or 
l0uetal.l [2^. Eor the sake of clarity we deleted the peaks 
of the spin period and its harmonics in Eig.|3 Here we are in¬ 
terested in the aperiodic variability of V0332-I-53 only, that is, 
in the evolution of the broad-band noise and QPOs as the out¬ 
burst decays. The coherence pulsation and its harmonics where 
fitted with narrow Lorentzians. The fundamental and the first 
two harmonics normally occupy one frequency bin. The third 
and fourth harmonic may not be resolved in some of the power 
spectra. 

In addition to the pulse noise, the power spectra of 
V0332-I-53 contains band-limited noise (that may turn into 
peaked noise) and QPO noise. The band-limited noise is repre¬ 
sented by three broad components that we shall call very-low 
frequency noise (Lvlfn). low-frequency noise (Llfn) and high- 
frequency noise (Lhfn)- The QPO noise consists of two com¬ 
ponents at 0.05 Hz (Llfqpo) and 0.22 Hz (Lhfqpo)- Only four 
Lorentzians (excluding the pulse noise) are present at any given 


time, though. Lvlfn is a zero-centred Lorentzian and describes 
the noise below ~ 0.02 Hz. Lvlfn and Lhfn are present both, 
in the high/soft state (HSS) and in the low/hard state (LHS), 
while Llfn is exclusive of the LHS. In the LHS the QPOs are 
absent. 

Table|2lgives the results for the power spectral fitting pro¬ 
cedure. 

3.3.1. Noise in the high/soft state (HSS) 

The soft branch comprises the intervals A-B-C-D. Interval E 
can be considered as a transition state between the HSS and the 
LHS. It would correspond to the so-called apex in Z sources. 
As the flux decreases the nns of Lvlfn increases, (from ~ 5% 
in region A to 12% in region E), while the nns of Lhfn shows 
an overall decreasing trend. QPO noise is present throughout 
the decay of the outburst but reaches its maximum significance 
at intermediate fluxes. A QPO, i.e., peaked noise with a quality 
factor Q - vjFWHM > 0.5, centred at ~ 0.2 Hz is weakly 
present in regions A and B and reaches maximum strength in 
region C and D. A second QPO at ~ 0.05 Hz emerges in re¬ 
gion C (it is only marginally present in regions A and B) and 
increases its power as the flux decreases, reaching maximum 
strength in region E with a fractional nns of 13%. The central 
frequencies of these QPOs do not change significantly. Note 
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Table 2. Power spectral parameters. Errors are 90% confidence level and upper limits 95%. 


Region 

Tvlfn 

Thfn 

Llfqpo 

Thfqpo 

Tlfn 

^Vdof 

Central frequency (Hz) 

A 

0 

0 

^•^^-0.08 

-0.05 

-0.03 


1.2/24 

B 

0 

fj 1 o-fO.05 

-0.05 

^•^'^-0.02 


1.7/24 

C 

0 

-0.01 

0-048!°;«“ 

0 991 +0.003 

u.z.z,i_ooo4 


1.2/21 

D 

0 

0.34t«| 

0-05 

0 999+0.003 
'-'•^^^-0.003 


2.0/24 

E 

0 

0.38t“ 

0-047!°:“‘ 


0 09+®-^^ 
'-'•^^-0.02 

2.3/20 

F 

0 

0-22tA? 



'-'•^^-0.02 

0.8/26 

Quality factor (vlFWHM) 

A 


0.1 


2.1 



B 


0.2 


4.6 



C 


0.4 

0.8 

1.2 



D 


0.4 

1.0 

1.8 



E 


0.5 

1.2 


0.3 


F 


0.8 



0.7 


Fractional rms 

A 

5.7 ±0.9 

13.3 ±0.7 

< 1.6 

2-8!^:) 



B 

6.1 ±0.6 

10 6+*^ ^ 
iu.u_07 

< 1.5 

2-6!i:I 



C 

7.2 ± 0.7 

8.8 ±0.1 

5.3 ± 0.3 

7 Q+0-4 
'•^-0.2 



D 

9.0 ± 1.1 

9.7 ±0.1 

6.0 ±0.1 

7.9 ± 0.2 



E 

12.0 ± 1.2 

6.2 ± 0.2 

12.8 ±0.3 


16.6 ±0.6 


F 

50 ±6 

24)^ 



40 ±6 



also, that the best-fit spectral parameters of the 0.22 HZ QPO 
given in Table |2 might be affected by the coupling of the QPO 
and the spin period components. 

3.3.2. Noise in the low/hard state (LHS) 

The study of the aperiodic variability in the hard branch is 
hampered by low statistics. The hard branch correspond to the 
lower count rates and it is reached when the 3-40 keV flux goes 
below ~ 5 X 10“® erg cm“^ s“' {Lx ~ 3 x 10^’ erg s“'). The 
overall 0.01-1 Hz fractional rms amplitude in the LHS is con¬ 
siderably larger than in the HSS. rms in region F is ~ 50%, 
compared to, for example, region C where rms ~ 10%. 

The 0.22 Hz QPO (Lhfqpo) that characterised the HSS dis¬ 
appears in the LHS. The 0.05 Hz QPO (Llfqpo) is still present 
in region E but not in region F. Instead, a new peaked noise 
component shows up, namely Llfn, which occupies the place 
left by the QPO noise. It first appears in interval E, but be¬ 
comes very strong in interval F. Lhfn, which during the HSS 
had Q < 0.5, becomes narrower in the LHS, with Q > 0.5. The 
pulse noise is highly suppressed, which may be a consequence 
of the overall increase in the strength of the continuum. 


3.4. Phase-lag spectra 

Since the power spectra has the disadvantage of ignoring phase 
information, we also calculated the cross spectrum between 
two concurrent light curves with different energies. If s{t) and 
h{t) are two such light curves and S (v) and H{v) denote their 
Fourier transform, then the power spectra of these light curves 
are respectively S *{v)S (v) and H*{v)H{v), while the cross spec¬ 
trum is defined as C(v) = S*{v)H{v). The * indicates the com¬ 
plex conjugate. The Fourier phase lag is the phase, i.e. the argu- 


£ 



Frequency (Hz) 


Fig. 4. Phase-lag spectra for each RXTE spectral region. The 
vertical dashed line marks the spin frequency (0.228 Hz). The 
lags correspond to the delay of 7-20 keV photons with respect 
to 2-7 keV photons. 


ment, of the complex quantity C(v). We obtained a mean phase- 
lag spectrum for each region of the CD (Fig.|3}. The lags were 
calculated between the 7-20 keV photons (hard band) with re¬ 
spect to 2-7 keV photons (soft band). 

Although the phase-lag behaviour of V0332-F53 as the out¬ 
bursts decays is rather complex, showing both, positive (hard 
photons lag soft photons) and negative lags (vice versa), soft 
lags dominate at high frequencies. Below ~0.1 Hz, no lags 
seem to be present. As in low-mass neutron-star and black- 
hole binaries, it is the phase lags, and not the time lags, that 
are roughly constant with Fourier frequency. 

Figure|3shows the average time lag as a function of the 4- 
30 keV flux for four frequency ranges. Two significant results 
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can be drawn; i) since the phase lags are constant, the aver¬ 
age amplitude of the time lags decreases with frequency and 
ii) at frequencies above the spin frequency the amplitude of the 
lags increases as the X-ray flux increases, i.e., as the source be¬ 
comes softer (the hard colour decreases). At lower frequencies 
the lags do not show any particular trend. 

4. Discussion 

The X-ray behaviour of transient Be/X-ray binaries is charac¬ 
terised by two types of outbursting activity; i) Type 1 outbursts 
are regular and (quasi)periodic, normally peaking at or close 
to periastron passage of the neutron star. The X-ray flux in¬ 
creases by about one order of magnitude with respect to the 
pre-outburst state, reaching L-c < 10^^ erg s“', ii) Type II 
are unexpected giant outbursts with increases in luminosity of 
~ 10^ - 10"^. They tend to last for several orbits and not neces¬ 
sarily coincide with periastron passage. 

Although aperiodic var iability studies of HMXB s have 
been carried out in the past (iBelloni & HasingerllQQflh . this is 
the first time that the aperiodic variability of a Be/X-ray binary 
has been systematically studied in correlation with the spectral 
evolution in the CD during a type II outburst. For a discussion 
on the origin of the QPOs and the variation of their spectral 
paramet ers with flux dur ing this same outburst the reader is re¬ 
ferred to lOu et ^ (l2005l) . 

4.1. Comparison with iow-mass X-ray binaries 

Two spectral states can be distinguished in the colour-colour di¬ 
agram (CD); a hard branch that corresponds to a low-intensity 
state (LHS) and a soft branch that corresponds to a high- 
intensity state (HSS). The hard branch is characterised by 
a low-amplitude change of the hard colour compared to the 
change in the soft colour. In the soft branch the amplitude of 
variability of the hard colour is about three times larger than 
that of the soft colour. The result is a pattern in the CD that 
resembles that of the Z sources. Like in Z sources the source 
moves gradually in the CD without jumps. However, there 
are also some differences which mainly affect the typical time 
scales and X-ray flux variations. Z sources trace out the Z track 
on time scales of hours to a day, while the motion of V0332-I-53 
in the CD is of the order of (tens of) days. Also, the amplitude 
of the X-ray lumin osity change over the Z trac k is typically less 
than a factor of 2 jdi Salvo et aiJl2nnnH2nn2h while the X-ray 
luminosity changes throughout the CD of V0332-I-53 are about 
a factor 60. 

In addition to the noise due to the X-ray pulsations we 
have identified three broad noise components and two QPOs in 
the 1/256^ Hz power spectra of V0332-I-53. The QPO noise 
is somewhat reminiscent of the horizontal (HBO) and normal 
branch (NBO) oscillations in Z sources; the 0.05 Hz QPO, be¬ 
ing only present near the hard branch or LHS (analogue to the 
’’horizontal branch”), would correspond to the HBO. The 0.2 
Hz QPO, whose frequency is independent of the position in 
the soft branch or HSS (analogue to the ’’normal branch”) and 
it is strongest in the middle of it, would be the counterpart to 
the NBO. Also, as in Z sources, the LFN is very weak or ab¬ 


sent in the HSS (’’normal branch”) but very strong in the LHS 
(’’horizontal branch”). Note, however, that although we have 
designated these components with the traditional names used 
in Z sources, namely, very-low (VLFN), low (LFN) and high- 
frequency noise (HFN), the time scales involved in V0332H-53 
are much longer than in Z sources. In particular, the central 
frequency of the LFQPO is three orders of magnitudes lower 
than typical values of HBOs, and the central frequency of the 
HFQPO is one order of magnitude lower than NBO values. 
Other difference between V0332-I-53 and Z sources is the vari¬ 
ation of the nns amplitude of the VLFN. In Z sources it de¬ 
creases as the source approaches the horizontal branch, while 
in V0332-I-53 the opposite is seen. 

_ h has been found jMiino et alJ l2002t loierlinski & Donj 

l2002h that low-mass neutron-star atoll sources that display 
large amplitude intensity variations {F„jaxlFmm~i00) also trace 
out three-branch patterns in the colour-colour diagram simi¬ 
lar to those of Z sources. However, the analysis of the rapid 
aperiodic variability in relation to the different regions of the 
colour-colour diagram and the actual time scales of the motion 
of the source through the diagram have re vealed very differ¬ 
ent behaviour between atoll and Z sources (Ivan Straaten et alJ 
l2003HReig et a iJ2nn4 . Here we And a similar result in a high- 
mass X-ray binary, namely, the source shows a Z track in the 
CD but the details of the motion in this track and time scales 
involved differ from those of low-mass Z sources. It seems that 
the Z topology might be a signature of the presence of a neu¬ 
tron star (hard surface) and that the details of the behaviour 
along this track depend on other factors like the type of opti¬ 
cal companion, strength of the magnetic held or mass transfer 
mechanism and deposition. 

Due to the transient nature of V0332-(-53, the unpredictabil¬ 
ity of type II outbursts and operational constraints, the observa¬ 
tions (both RXTE and INTEGRAL) covered the decay of the 
outburst only. It would be very interesting to study the transi¬ 
tion from low-intensity to high-intensity states and whether the 
source follow the reverse path in the CD or there is some degree 
of hysteresis (i.e., the fact that the transition between states in 
one direction does not occur for the same value of the spectral 
and timing parameters as in the rev erse direction) as have been 
see n in low-mass atoll s ources (e.g. lReig et al .12003) and BHC 
(e.g iBelloni et alj|2005l) . Note that hysteresis is not present in 
Z sources. 

4.2. Coupiing of the periodic and aperiodic noise 
components 

iLazzati & Stellallll997l) investigated the relationship between 
the periodic and aperiodic variability in three accreting X-ray 
pulsars and found that a highly significant coupling is present. 
These authors used a shot noise model to account for the aperi¬ 
odic component and argued that if the shots represent inhomo¬ 
geneities in the accretion flow (blobs or clumps of matter) and 
are produced close to the neutron star surface then they should 
be also modulated by the same mechanism that gives rise to the 
periodic signal. This coupling manifests itself by the broaden¬ 
ing of the wings of the narrow peaks due to the periodic mod- 
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Fig. 5. Time lags as a function of the X-ray flux for four fre¬ 
quency ranges. 



Frequency (FIz) 


Fig. 6. Coupling of the periodic and aperiodic variability, evi¬ 
denced by the broadening of the wings of the spin period peak. 
The power of RXTE regions C and E was shifted for plotting 
purposes. 


ulation and becomes more apparent when the red-noise power 
increases short-ward of the pulsar frequency. The width of the 
wings are strongly dependent of the ratio of the pulsar period P 
to the characteristic decay time of the shots t. If 2jiP > r then 
the wings are so broad that they become indistinguishable from 
the power spectra continuum. 

The 0.22 Hz QPO detected in V0332-F53 constitutes the 
first detectio n of a QPO ri ding on the spin frequency of a 
neutron star don et al .ll2nnli and points toward a strong cou¬ 
pling between the periodic and red-noise components. Pigure|3 
shows the power spectra of V0332-F53 corresponding to var¬ 
ious spectral states without the removal of the peaks that re¬ 
sult from the spin period and its harmonics. The coupling be¬ 
tween the periodic and aperiodic variability is strong during 
intervals C to E and weak in intervals A and B. Since the spin 
period does not change substantially, the characteristic decay 
time of the shots must be significantly smaller at the beginning 
of the outburst decay. If the shots arise from blobs in the accre¬ 
tion flow, and if the shot life time is associated with the blob 


size then we conclude that at high luminosities the blobs are 
smaller than at low luminosities. Maybe the radiation pressure 
breaks the accretion flow into small discrete clumps of matter. 
Alternatively, the weak coupling in intervals A and B could be 
due to the flatter contin uum short-ward of the spin frequency 
(iLazzati & StellJl997h . 


4.3. Phase lags 


Eigure shows the time lags as a function of X-ray flux for 
four frequency ranges. If we assume the interpretation that each 
Eourier component is emitted by a specific spatial scale in the 
accretion flow then the highest frequencies would correspond 
to the closest distances to the compact object. The different be¬ 
haviour of the lags depending on the frequency range consid¬ 
ered allows us to put some constraints on the size of the accre¬ 
tion column. At high X-ray flux the amplitude of the time lags 
in the frequency range 1-8 Hz is ~ 3 ms, which measured as 
light crossing time corresponds to a length of scale of ~ 9x10^ 
cm or ~ 0.3 magnetospheric radii. We have assum ed the radius 
of the magnetosphere to be 2.9 x 10^ cm (see e.g. lFrank et alJ 
Il992il . As the flux decreases, i.e., as the source spectrum hard¬ 
ens, the magnitude of the lags decreases, which would imply 
that the X-ray emitting region (presumably the accretion col¬ 
umn) reduces its size. At the end of the outburst, close to quies¬ 
cence, the lag is consistent with zero, which can be interpreted 
as the disappearance of the accretion column. An increase of 
the magnitude of the lags (associated with QPO noise) as the 
flux increases ha s been also repo rted for the black-hole system 
XTE J1550-564 jCiii et alJboncIl . Further support in favour of 
the reduction in size of the accretion column as the source gets 
fainter is provided by the fact that the c yclotron resonance en - 
ergy increases during the burst decline (|^towlavi^e^ayj200^. 
This is in agreement with Mihara’s model dMihMTeTalEoO^) 
that predicts this behaviour as due to the reduction of the ac¬ 
cretion column scale height as a function of source intensity. 


At lower frequencies there is not a clear pattern of vari¬ 
ability of the lags versus the X-ray flux. In fact, below the 
frequency that corresponds to the pulse peak, no significant 
time delays are present in the data for most frequencies. Low 
frequencies sample photons produced further away from the 
compact object, perhaps in the region where the accretion flow 
interacts with the magnetosphere. A time lag of 0.05 s corre¬ 
sponds to a length scale of ~ 1.5 x 10® cm or ~ 5 magneto- 
spheric radii. The lack of significant lags at low frequencies 
may also indicate that most of the X-ray radiation is produced 
in the vicinity of the neutron star. 


The phase- l ag spe ctra of neutron-star (iFord et al.l Il999t 
|01iYe,&^ane]J_ 200 ih and Galactic b lack-holes binaries 
dNowak et alJll999Hwiinands et al.lll999li are very similar but 
incompatible with the assumption of a lag constant in time as 
predicted by Comptonisation in a uniform medium. Instead, 
observations show the phase lag to be roughly independent 
of Fourier frequency. The same result is now seen in the ac¬ 
creting X-ray pulsar V0332-F53. In V0332-F53, however, soft 
lags dominate the 0.1-10 Hz phase-lag spectra. Soft lags are 
not a rare phenomenon. They have been observed in neutron- 
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star 1 

Ou et al.ll2004l) and black-hole systems JWiinands et alJ 

|l999 

Cui et al.l2000HLin et alJl200(ilReig et aljboooll. While 


models based on Comptonisation of low-energy photons in a 
non-uniform corona can reproduce the amplitude and energy 
dependence of the hard lags the change in sign of the lags 
are hard to explain by invoking light-travel ti me differences 
betwe en photon s at different en ergies (but see lOhkawa et alJ 
I 2 OO 5 I) . Recently. IVarnier j ( l2005l) has proposed that partial ab¬ 
sorption (related to the jet mechanism) of the signal can ac¬ 
count for the phase lag reversal in microquasars. Since X-ray 
pulsars do not show radio jets, other mechanisms producing the 
absorption of the soft/low-energy X-rays must be at work if the 
phase lag behaviour of V0332-H53 is to be explained with this 
model. 

5. Conclusion 

The study of the evolution of the X-ray colours of the high- 
mass X-ray pulsar V0332H-53 through the decay of a type II 
outburst has shown that these type of systems can also trace 
out Z-shaped pattern in the colour-colour diagram, like its low- 
mass cousins, making this feature a possible signature of the 
presence of a neutron star. While colour-colour diagrams may 
represent a useful tool for the detection of spectral states, they 
cannot unveil the true nature of the source. Analysis of the rapid 
aperiodic variability, allowing the study of the noise compo¬ 
nents in different regions of the colour-colour diagram and the 
actual time scales of the motion of the source through the dia¬ 
gram, is crucial to identify source type and state. Unlike black- 
hole binaries and atoll sources but similarly to the Z source 
GX 5-1, the phase-lag spectra of the noise below 10 Hz in 
V0332H-53 are dominated by soft lags. We conclude then that 
although the mechanism responsible for the lags does not de¬ 
pend on the presence or absence of a hard surface it may be 
affected by the magnetic field, which in turn may dictate the 
geometry of the accretion flow, that is, whether an accretion 
column is formed. 
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